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1. Introduction 

It is expected that strongly interacting matter undergoes a transition in some temperature in- 
terval from hadron gas to deconfined state also called the quark gluon plasma (QGP) [|l|]. Creating 
deconfined medium in a laboratory is the subject of the large experimental program at RHIC [§] 
and is going to be the goal of the future heavy-ion program at LHC Attempts to study QCD 
thermodynamics on the lattice go back to the early 80's when lattice calculations in SU (2) gauge 
theory provided the first rigorous theoretical evidence for deconfinement ||, The problem 
of calculating thermodynamic observables in pure gluonic theory was solved in 1996 while 
calculation involving dynamical quarks were limited to large quark masses and had no control over 
discretization errors |j8|, ^, [T^ (see Refs. [|n|, 12] for reviews). Lattice calculations of QCD ther- 



modynamics with light dynamical quarks remained challenging until recently. During the past 5 
years calculations with light u,d quarks have been performed using improved staggered fermion 
l3 111, ll, O, [11, llll 1^ III, II g, II ||] (see Refs. g.) 



actions 



To get reliable predictions from lattice QCD the lattice spacing a should be sufficiently small 
relative to the typical QCD scale, i.e. Aqcd<^ ^ 1- For staggered fermions, which are used for 
calculations at non-zero temperature, discretization errors go like ^((aAgcz))^) but discretization 
errors due to flavor symmetry breaking turn out to be numerically quite large. To reduce these 
errors one has to use improved staggered fermion actions with so-called fat links []30|]. At high 
temperature the dominant discretization errors go like {aT)-^ and therefore could be very large. 
Thus it is mandatory to use improved discretization schemes, which improve the quark dispersion 
relation and eliminate these discretization errors. Lattice fermion actions used in numerical calcu- 
lations typically implement some version of fat links as well as improvement of quark dispersion 
relation and aie referred to as p4, asqtad, HISQ and stout. In lattice calculations the temperature is 
varied by varying the lattice spacing at fixed value of the temporal extent A^^- The temperature T is 
related to lattice spacing and temporal extent, T = 1/ {Nxo). Therefore taking the continuum limit 
corresponds to A^^ — ^ °° at the fixed physical volume. For the same reason discretization errors in 
the hadronic phase could be large when the temperature T is small. 

In this paper I discuss lattice QCD calculation of the transition temperature, the equation of 
state as well as different spatial and temporal correlation functions. 



2. Equation of State 

The equation of state has been calculated with pA and asqtad action on lattices with temporal 
extent Nx = A, 6 and 8 |l8|, |^]. In these calculations the strange quark mass was fixed to 
its physical value, while the light (m, d) quark masses 10 times smaller than the strange quark 
mass have been used. These correspond to pion masses of (220 — 260) MeV. The calculation of 
thermodynamic observables proceeds through the calculation of the trace of the energy momentum 
tensor e — 3p also known as trace anomaly or interaction measure. This is due to the fact that this 
quantity can be expressed in terms of expectation values of local gluonic and fermionic operators. 
The explicit expression for £ — 3;? in terms of these operators for p4 and asqtad actions can be 



found in Ref. [20|. Different thermodynamic observables can be obtained from the interaction 
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measure through integration. The pressure can be written as 

p{T) p{To) f^dT' 



7^4 



To' 



(2.1) 



The lower integration limit To is chosen such that the pressure is exponentially small there. Fur- 
thermore, the entropy density can be written as s = {e + p)/T. Since the interaction measure is the 
basic thermodynamic observable in the lattice calculations it is worth discussing its properties more 
in detail. In Fig. |I| I show the interaction measure for p4 and asqtad actions for two different lat- 
tice spacings corresponding to A^^ = 6 and 8. In the high temperature region, T > 250 MeV results 
obtained with two different lattice spacings and two different actions agree quite well with each 
other. Furthermore, recent calculations with HISQ actions also give results for the trace anomaly 



which are consistent with these [|24|]. Discretization errors are visible in the temperature region, 
where e — 3p is close to its maximum as well as in the low temperature region. At low tempera- 
tures the lattice data have been compared with the hadron resonance gas (HRG). As one can see 
the lattice data fall below the resonance gas value. This is partly due to the fact that the light quark 
masses are still about two times larger than the physical value as well as to discretization errors. 
The ^((qAqcd)^) discretization errors in the hadron spectrum are suppressed at high temperatures 
as the lattice spacing a is small there. Also hadrons are not the relevant degrees of freedom in this 
temperature region. But at small temperatures, where hadrons are the relevant degrees of freedom, 
these discretization effects are significant. It turns out, however, that the HRG model that takes into 
account the quark mass dependence and discretization errors in the hadron spectrum can describe 



the lattice data quite well [|31|]. The large discretization errors in the low temperature region is the 
reason for the discrepancies with the stout results at temperatures T < 200 MeV [^J]. The lattice 
results for the trace anomaly obtained with stout action are also different in the high temperature 
limit 1 23]. Namely, {e — 3p)/T'^ calculatied with stout action is about 50% smaller in the high 
temperature region. To resolve this problem calculation at finer lattice spacing are needed. 
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Figure 1: The interaction measure calculated with p4 and asqtad actions in entire temperature range (left) 
and at low temperatures (right) from Ref. [pO[]. The dashed and dashed-dotted lines are the prediction of 
hadron resonance gas (HRG) with all resonances included up 2.5GeV (dashed) and 1 .5GeV (dashed-dotted), 
respectively. 

The pressure, the energy density and the entropy density are shown in Fig. ^. The energy 
density shows a rapid rise in the temperature region (185 — 195) MeV and quickly approaches 
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about 90% of the ideal gas value. The pressure rises less rapidly but at the highest temperature it 
is also only about 15% below the ideal gas value. In the previous calculations with the p4 action it 
was found that the pressure and energy density are below the ideal gas value by about 25% at high 



temperatures [|10p. Possible reason for this larger deviation could be the fact that the quark masses 
used in this calculation were fixed in units of temperature instead being tuned to give constant 



meson masses as lattice spacing is decreased. As discussed in Ref. [ pq ] this could reduce the 
pressure by 10 — 15% at high temperatures. In Fig. ^I also show the entropy density divided by 
the corresponding ideal gas value and compare the results of lattice calculations with resummed 
perturbative calculation [32, ^ as well as with the predictions from AdS/CFT correspondence for 



the strongly coupled regime [34]. The later is considerably below the lattice results. Note that 
pressure, energy density and the trace anomaly have also been recently discussed in the framework 



of resummed perturbative calculations which seem to agree with lattice data quite well [|35|]. 

The differences between the stout action and the p4 and asqtad actions for the trace anomaly 
translates into the differences in the pressure and the energy density. In particular, the energy 
density is about 20% below the ideal gas limit for the stout action. 
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Figure 2: The energy density and the pressure as function of the temperature (left), and the entropy density 
divided by the corresponding ideal gas value (right). The dashed lines in the right panel correspond to the 
resummed perturbative calculations while the solid black line is the AdS/CFT result. 



3. Chiral and deconflnement transition 

The finite temperature transition in QCD has aspects related to deconfinement and chiral sym- 
metry restoration. Deconfinement aspects of the transition are related to color screening and, for 
infinitely heavy quarks, also to the center symmetry. The order parameter for deconfinement is 
the Polyakov loop which is related to the free energy of isolated static quark. The breaking of 
the center symmetry is signaled by non-zero value of the Polyakov loop. In the opposite limit of 
massless quarks QCD has the chiral symmetry. The quark condensate {'^f'^f) is the order parameter 
for this symmetry. The chiral symmetry is broken in the vacuum and expected to be restored at 
high temperatures. In Fig. ^ the renormalized Polyakov loop and the subtracted chiral condensate 
are shown. Both quantities show a smooth change in the transition region which is consistent with 
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Figure 3: The subtracted chiral condensate (left) and the renormalized Polyakov loop (right) calculated for 
different actions. In the left panel we take Tc — 164 MeV for the chiral temperature in QCD. The numerical 
data are taken from ^ g ||, |^ . 



the fact that the finite temperature transition is an analytic crossover and not a true phase transi- 



tion []37[]. The universal aspects of the chiral transition seem to be relevant for the range of quark 
masses which include the physical light quark mass and the corresponding transition temperature 
can be defined [p8|]. Calculations with the p4 action on A^^ = 4 and 6 lattices gave an estimate 
Tc = 192(4) (7)MeV for the chiral transition temperature [15] wliich is significantly larger than the 
value of about 155MeV obtained for stout action ^ This is due to large flavor symmetry break- 
ing for the p4 action. As one can see from Fig. ^ the discrepancies between different actions 
are reduced when considering larger A^^ or the HISQ action where the effects of flavor symmetry 
breaking are much smaller. The preliminary estimate based on asqtad action for the chiral transition 
temperature based on asqtad action is Tc = (164 it 6) MeV [^9|]. The effects of flavor symmetry 
violations are smaller for the renormaUzed Polyakov loop. Only for temperatures T < 180 MeV 
there are some discrepancies between the p4 and asqtad results obtained on A^^- = 8 lattices and 
the stout results. The asqtad results on A^^. = 12 lattices as well the HISQ results are in very good 
agreement [24, ^]. In Fig. || the renomalized Polyakov loop is shown for HISQ and stout actions 
as function of T /Tc. The results are compared to the renormalized Polyakov loop in pure gauge 
theory. Here Tc denotes the phase transition temparature for pure gauge theory and the chiral tran- 
sition temperature for QCD. As one can see except for very high temperatures the Polyakov loop 
in QCD is very different from the one in pure gauge theory and therefore it is not clear if the center 
symmetry plays an important role in QCD with light dynamical quarks. 



4. Correlation functions of static quark anti- quark pair 

One of the most prominent feature of the quark gluon plasma is the presence of chromoelectric 
(Debye) screening. The easiest way to study chromoelectric screening is to calculate the singlet 

' Due to the crossover nature of the transition different chiral observable results in different estimate of the transition 
temperature, the so-called renormalized susceptibility calculated with stout action gives Tc ~ 147 MeV 
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Figure 4: The singlet free energy Fi{r,T) calculated in Coulomb gauge on 16^ x 4 lattices (left) and the 
combination Fi{r,T) ~ Foo{T) as function of rT (right). The solid black line is the parametrization of the 
zero temperature potential. 



free energy of static quark anti-quark pair (for recent reviews on this see Ref. pT| , ^]), which is 
expressed in terms of correlation function of temporal Wilson lines in Coulomb gauge 



exp(-Fi(r,r)/r) = lTr(W(r)W^(0)). 



(4.1) 



L = TrW is the Polyakov loop. The singlet free energy is in fact the logarithm of the static meson 
correlation function evaluated at T = l/T. Instead of using the Coulomb gauge the singlet free 
energy can be defined in gauge invariant manner by inserting a spatial gauge connection between 
the two Wilson lines. Using such definition the singlet free energy has been calculated in SU (2) 
gauge theory [p3|]. It has been found that the singlet free energy calculated this way is close to the 
result obtained in Coulomb gauge [^]. The singlet free energy turned out to be useful to study 
quarkonia binding at high temperatures in potential models (see e.g. Refs. [ p3| , ^ ^ ^ p9|]). 
The singlet free energy also appears naturally in the perturbative calculations of the Polyakov loop 
correlators at short distances [p]]. 

The singlet free energy was recently calculated in QCD with one strange quark and two light 
quarks with masses corresponding to pion mass of 220MeV on 16^ x 4 lattices [|50|]. The numerical 
results are shown in Fig. Q At short distances the singlet free energy is temperature independent 
and coincides with the zero temperature potential. In purely gluonic theory the free energy grows 
linearly with the separation between the heavy quark and anti-quark in the confined phase. In 
presence of dynamical quarks the free energy is saturated at some finite value at distances of about 
1 fm due to string breaking 53 1. This is also seen in Fig. ^ Above the deconfinement 



temperature the singlet free energy is exponentially screened at sufficiently large distances [ |40| , |52| ] 
with the screening mass proportional to the temperature , i.e. 



Fi(r,r)=fec(r) 



exp(— m/)(r)r), ~ T. 



(4.2) 



3 47ir 

Therefore in Fig. ^ we also show the combination Fi (r, T) — Foo{T) as a function of rT. As one can 
see from the figure this function shows an exponential fall-off at distances rT > 0.8. The fact that 
the slope is the same for all temperatures means that mjj ~ T, as expected. 
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5. Heavy and light meson correlation functions 

Information on hadron properties at finite temperature as well as the transport coefficients are 
encoded in different spectral functions. In particular the fate of different quarkonium states in QGP 
can studied by calculating the corresponding quarkonium spectral functions (for a recent review see 



Ref. [42]). On the lattice we can calculate correlation function in Euclidean time. This is related 



to the spectral function via integral relation 

G(T,r) = rdcoa{co,T)K{T,co,T), K{t,co,T) = • (5-1) 

Jo smh(G)/2y j 

Given the data on the Euclidean meson correlator G{r,T) the meson spectral function can be 
calculated using the Maximum Entropy Method (MEM) [^]. For charmonium this was done by 



using correlators calculated on isotropic lattices [55, pq] as well as anisotropic lattices [57, 58 



59] in the quenched approximation. It has been found that quarkonium correlation function in 



Euclidean time show only very small temperature dependence [56, 59]. In other channels, namely 
the vector, scalar and axial- vector channels stronger temperature dependence was found [ p6[ ^ ]. 
The spectral functions in the pseudo-scalar and vector channels reconstructed from MEM show 



peak structures which may be interpreted as a ground state peak [^, ^ Together with the 
weak temperature dependence of the correlation functions this was taken as strong indication that 
the IS charmonia (T]^ and J/y) survive in the deconfined phase to temperatures as high as 1.67^ 



[ ]56[ , |57| , |58| ]. A detailed study of the systematic effects show, however, that the reconstruction of the 
charmonium spectral function is not reliable at high temperatures [|^], in particular the presence 
of peaks corresponding to bound states cannot be reliably established. Presence of large cutoff 
effects at high frequencies also complicates the analysis [|62[]. The only statement that can be made 
is that the spectral function does not show significant changes within the errors of the calculations. 
Recently quarkonium spectral functions have been studied using potential models and lattice data 
for the singlet free energy of static quark anti-quark pair [ p7| , 48, ^]. These calculations show 



that all charmonium states are dissolved at temperatures smaller than 1.2rc, but the Euclidean 
correlators do not show significant changes and are in fairly good agreement with available lattice 
data both for charmonium [ ^ , ^ and bottomonium [ ]59| , |^]. This is due to the fact that even 
in absence of bound states quarkonium spectral functions show significant enhancement in the 
threshold region [p6|]. Therefore previous statements about quarkonia survival at high temperatures 
have to be revisited. Exploratory calculations of the charmonium correlators and spectral functions 



in 2-fiavor QCD have been reported in Ref. [ ]61[ ] and the qualitative behavior of the correlation 
functions was found to be similar. 

The large enhancement of the quarkonium correlators above deconfinement in the scalar and 



axial-vector channel can be understood in terms of the zero mode contribution [ [46[ , |63|] and not 
due to the dissolution of the IP states as previously thought. Similar, though smaller in magni- 
tude, enhancement of quarkonium correlators due to zero mode is seen also in the vector channel 



[p9[]. Here it is related to heavy quark transport [^, Due to the heavy quark mass the Eu- 
clidean correlators for heavy quarkonium can be decomposed into a high and low energy part 
G{r,T) = Giow(i',r) + Ghigh(T,r) The area under the peak in the spectral functions at zero en- 
ergy ft) ~ giving the zero mode contribution to the Euclidean correlator is proportional to some 
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susceptibility, G\^j^^{r,T) ~ Tx'{T), which have been calculated on the lattice in Ref. []65[]. It is 
natural to ask whether the generalized susceptibilities can be described by a quasi-particle model. 
The generalized susceptibilities have been calculated in Ref. [ ]66| ] in the free theory. Replacing the 
bare quark mass entering in the expression of the generalized susceptibilities by an effective tem- 
perature dependent masses one can describe the zero mode contribution very well in all channels 

While temporal correlators are not sensitive to the change in the spectral functions spatial 
quarkonium correlation functions could be more sensitive to this. Recent lattice calculations show 
indication for significant change in spatial charmonium correlators above deconfinement [|67|] . 

The spectral function for light mesons has been calculated on the lattice in quenched approx- 
imation [ ]68| , |69[ , [70[ , [7T| ]. However, unlike in the quarkonia case the systematic errors in these 
calculations are not well understood. The spatial meson correlation function provide interesting 
insight into the modification of the spectral functions also in the light quark sector. In particular, 
the study of the spatial meson correlation function in the vector and axial vector channel indicate 
the degeneracy of the vector and axial vector spectral functions at the chiral transition |72]. At 
the same time the lattice calculations of the pseudo-scalar and scalar correlation function indicate 
that the ?7a ( 1 ) axial symmetry is only restored at temperatures significantly higher than the chiral 
transition temperature [72|. 



6. Conclusions 



In recent years significant progress has been achieved in studying strongly interacting matter 
at high temperatures using lattice QCD. Equation of state and transition temperature have been 
calculated at several lattice spacings allowing for controlled continuum extrapolations. Current 
lattice data suggest a transition temperature for chiral symmetry restoration of (147 — 164) MeV. 
Temporal meson correlation functions have been also studied in lattice QCD but the numerical 
data are not precise enough to provide detailed information of the corresponding meson spectral 
functions. However, these lattice data are useful to constrain the model calculations of the meson 
spectral functions. Spatial meson correlators provide additional information about meson spectral 
functions, in particular indication of effective restoration of the [/a(1) symmetry at temperatures 
larger than the chiral transition temperature. 
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